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With the purpose of increasing the information about the lower frequency shielding 
effectiveness of closed structures composed of various arbitrary materials, an idealized 
problem is considered of the scattering and absorption of a plane electromagnetic wave 
impinging on a spherical shell. The region between the outer and inner radius of the sphere 
contains an arbitrary material; the other two regions are free space. The absorption and 
scattering coefficients for the incident plane wave are expressed as an infinite sum of spherical 
waves, using spherical Bessel functions. The coefficients are calculated numerically and 
shown for a large range of complex material parameters and frequencies from 100 kc/s to 
1000 Mc/s. Equations were programmed in FORTRAN, and numerical calculations ac- 
complished on the IBM 1620 computer. Combinations of complex material properties of 
permeability, permittivity, and conductivity were sought which gave large, broadband 
absorption coefficients. Several promising compositions have resulted, and are presented. 
The effects of varying the individual parameters are also discussed. 

1. Introduction 

In order to utilize a mathematically complete and general solution to the boundary value 
problem, the symmetrical configuration of a spherical shell has been used to approximate the 
closed walls of the normally cubical or rectangular shielded enclosure. The behavior of the 
impinging plane wave electromagnetic field external to the spherical structure is considered. 
In order to completely describe the perturbing effect of the spherical shell, knowledge of the 
field outside, inside, and within the shell material is necessary. This investigation has been 
primarily directed toward a systematic exploratory search for wall (or shell, in this case) 
materials that will give uniform absorption of plane EM waves in the broad spectrum between 
0.1 and 1000 Mc/s. 

The assumptions made in this investigation are: 

1 . The spherical structure is isolated in space, 

2. the wave incident on the surface of the outer sphere is a plane wave, 

3. the shell material is homogeneous and isotropic, 

4. the shell is composed of an arbitrary material with complex properties of /x*, e*, and a. 
The inner sphere and the exterior region is free space. 

The geometry of the spherical configuration is shown in figure 1. For this case, the equations 
describing the absorption (C a ) and scattering (C s ) coefficients have been derived (appendix 
A). The Mie scattering equations [Stratton, 1941] expanded to two concentric spheres by 
Salati [1962], are essentially those used in this paper. We have found that it is possible to 
extend the equations to cover three concentric spheres but only the two sphere solution is 
presented here. A general solution for any number of concentric spheres has been given by 
Wait [1962 and 1964]. 

The shell material parameters have been systematically varied over a wide range of values. 
Ranges of values of parameters were selected generally on the basis of known conductors or 
nonconductors. The case of high conductors yields considerable simplification of the C s and 
C a equations and has been treated separately. The dimensional properties of the spherical 
shell structure have, in most cases, been selected to approximate the dimensions of actual 
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Figure 1. Diagram of spherical configuration. 

Inner and outer regions are free space. 
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shielded enclosures. The majority of data was computed on the basis of a 2 (outer radius) = 5 
m (16.4 ft) and a shell thickness (a 2 — Gi) = 0.15 m (5.9 in.). In many cases data were taken at 
different radii to show the variation in C a and C s . The effect of shell thickness variation was 
also considered. 

Approximately 150 separate combinations of material parameters or dimensional variations 
were computed. Only a summary of the materials that gave large absolute values of absorption 
coefficients and broadband absorption is given here. 

2. Absorption and Scattering Coefficients for the Case of Arbitrary Materials 

The absorption coefficient, C a , is the ratio of the absolute absorption cross section to the 
geometrical cross section. It is also the ratio of the total power absorbed to the power of the 
incident wave crossing an area equal to the geometric cross section. Similar definitions hold 
for the scattering coefficient, O s . 

According to (A 15) and (A 16) in the appendix, the scattering and absorption coefficients 
may be expressed in terms of the external field coefficients, a T n and b r n , as 



c, 



=^S(2«+1)(K| 2 +|^| 2 ) 

£2 n-l 



0.= -SS(2«+l)(K| 2 +|6;| 2 +Reaj;+Ee6j). 

£2 11 = 1 

The field coefficients may be expressed in terms of first and second spherical Hankel func- 
tions, hn ] andA£ 2) and spherical Bessel functions of the first kmd,j n . The two external coeffici- 
ents have identical form except for the interchange of the relative admittance of the material 

(g=-y/K e /K m ) and the relative impedance (Z=- =-\JK m K^) . Writing only one of these 



coefficients, we have from (A7) 

[ 



jn(h) 
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where 

^=k a h £ 2 = k a 2 , pi = ka u p 2 =Jca, (k=k \K m K e ) 
and 

a< n (z)- A *~ l(a 



' Ai"(*) 


x' 


.7™_i(a") 
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j»(«) 


x' 


_A^i(«) 


n 



dn ^ = h%\x) X 

The material is characterized by complex permeability and permittivity, 

e*=e (K' e -iK' e ')=e Q K e 

where 

K , /=K f e loss tangent for dielectrics 

K' e ' = -- — for conductors. 

CO 

3. Dielectric Material 

In general, high dielectric materials will give absorption curves with values that are much 
less than unity and have severe fluctuations in magnitude over the frequency range. For high 
dielectric materials small changes in thickness or in the radius of the shell structure will cause 
drastic changes in the shape of the absorption curves. This effect is illustrated for a change 
in radius in figures 2 and 3, where a dielectric constant of 100 has been used. The radius of 
the spheres has been changed by a factor of two. The scattering coefficient for the same material 
is shown in figure 4. The numbers shown in parentheses are (K' e , K") (iC, K^), respectively. 

4. Magnetic Material 

The magnetic properties have been investigated on the basis of variation of the permeability 
(K' m ) between 1 and 100 (see fig. 5) and a magnetic loss factor (iC) between 0.1 and 100 (see fig. 
6). The above data is based on a dielectric loss factor {K' e ') of 100 and dielectric constant (K' e ) 
of one. Small changes in the general shape of the curves indicate that wide variations in K' m 
and K^ have little effect on the uniformity of C a . Figures 7 and 8 show the effect on C a for 
variations of conductivity between 10~ 4 and 10 4 mhos/m and a permeability of 100, using a 
loss tangent of one. The most uniform curve for C a occurs at (1.1) (100,100), indicating that a 
low conductivity (<r=l mho/m) is desirable in this combination. 

5. High Conductivity Materials 

The high conductivity case (<j>10 4 mhos/m) occurs whenever the skin depth is an order of 
magnitude smaller than the shell thickness, so that the fields are damped before reaching the 
inner sphere. For this case the field coefficients are, from (A14), 

, = j.(h) 4 1) U)-2y„fe ) 

734-668— 64 7 " ^7 























































OUTER RADIUS ■ 5 METER! 






























































WALL THICKNESS = 0.15 


































METERS 




































































h- 




















































UJ 

o 
fr io- 2 






















































« 






















































o 








































































































1 






































































































0. 




















































£ 






















































































































































































/ 




















































' 


















































/ 


















































/ 































FREQUENCY- Mc/S 



Figure 2. Absorption coefficient versus frequency in Mc/s 

for material with high dielectric constant (100,0.1) (1,0). 

a2 = 5m. 



10 -' 




pn 


i 


rim 


JS ■ 10 METERS 


i i 




















* 






= 


OUTER *ADl 


= 


























WALL THICKNESS = 015 METER 


S 
























































































t ,0 " 2 








































H 




















-^ 










I 




ii 
























1 
















g 




















f 




















^ 








































or 
o 




















1 




















< 
































































































/ 








































































| 








/ 




































/ 


<' 























FREQUENCY- MC/S 



Figure 3. Absorption coefficient versus frequency in Mc[s 
for material with high dielectric constant (100,0.1) (1,0). 
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Figure 4. Scattering coefficient versus frequency in Mc/s for 
material with high dielectric constant (100,0.1) (1,0). 
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FREQUENCY- MC/S 



Figure 5. Absorption coefficient versus frequency in Mc/s 
for variations in permeability between 1 and 100. 
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Figure 6. Absorption coefficient versus frequency in Mc/s 
for variations between 0.1 and 100 in the magnetic loss factor. 
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Figure 7. Absorption coefficient versus frequency in Mc/s, 

Variations of conductivity from 10-i to 10-* mho/m for high permeability ma- 
terial. ,u=100. 
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Figure 8. Absorption coefficient versus frequency in Mc/s. 

Variation of conductivity from 1 to 10* mho/m for high permeability material. 

M=100. 
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SCATTERING COEFFICIENT C s 

Figure 9. Scattering coefficient for conductive materials versus 
megacycle — meter product. 
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Using this, the scattering and absorption coefficients were computed for various material 
parameters. The resulting data has the following form: 

(1) The scattering coefficient is independent of the material parameters and is a function 
only of the product Fa 2 (see fig. 9). For Fa 2 less than 20 Mc/s-meter, the scattering coefficient 
varies as the fourth power of the product Fa 2 . This is Rayleigh scattering and corresponds to 
a wavelength greater than 10 times the radius. For Fa 2 greater than 50 Mc/s-meter, the 
scattering coefficient is approximately two. 

(2) The absorption coefficient varies as 



C 



{ ^ + f 2+K ™ J 2 UF,a 2) 



(2) 



where j a (F y a 2 ) is shown in figure 10 for radii of 1, 5, 10, and 20 m. This factor has only a 
slight dependence on the radius and may be approximated by 



/a^W^XlO- 3 . 

See figure 11 for data on high conductivity materials (o-=10 4 to o-=10 8 ). 



Figure 10. 
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Absorption factor versus frequency in Mc/s for 
conductive materials. 
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Figure 11. Absorption coefficient versus frequency for material 
conductivity variations between 10* and 10 l mho/m. 



Figure 12. Absorption coefficient versus frequency in Mc/s 
for material conductivity variations between 0.1 and 10~ A 
mho/m. 
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Figure 13. Absorption coefficient versus frequency in Mc/s 

for material of 0.1 mho/m conductivity. 

02= 5m, 20m. 
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Figure 14. Absorption coefficient for ferramic. 
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6. Conductors 

The approximations made in the high conductivity case are no longer valid for materials 
less than 10,000 mhos/m. For calculation of C a where <r<10 4 mhos/m it has been necessary 
to use the general equations describing any material. The results of low conductivities (1.0 
to 10 -4 mhos/m) are indicated in figure 12. The most uniform curve of any material examined 
is given by a material conductivity of 0.01 mhos/m. 

Variation of the absorption curve for a change in the radius of the sphere has been investi- 
gated. The effect is to shift the maxima downward in frequency as the radius is increased. 
Figure 13 shows the effect of a radius change from 5 to 20 m for material a- of 0.1 mhos/m. 
This effect is true for most low to medium conductivity materials. 

7. Variation of Material Parameters With Frequency 

Generally, the variation in K' e , K'J , K' m , K'J, for most materials is not great over the frequency 
range of interest in this study. However, the list of materials [Von Hippie, 1954] indicates 
there is a group of materials whose parameters vary rapidly in the range from 0.1 to 1000 
Mc/s. The listed parameter values for a ferramic were used to determine C a as shown in figure 
14. Even with the large changes in K e and K m the curve appears uniform. However, only 
the data points shown were available and there may be variations due to a lack of data in the 
intermediate ranges. 

8. Conclusions 

The uniform absorption curves occur for materials whose conductivity is in the region 0.1 
to 0.01 mhos/m. Some magnetic materials also give near-uniform absorption curves. Die- 
lectric materials with high dielectric constants and low loss factors are definitely not suitable. 
Low dielectric constant and high loss factor materials exhibit the characteristics of uniform 
absorption; however, the low conductivity materials make more suitable building materials. 

The absorption coefficient should be as independent as possible of variations in the dimen- 
sional properties of the structure. The choice of a material will depend not only on 
the uniformity of the absorption versus frequency curve, but also the behavior of the absorption 
coefficient with changes in radius and shell thickness. 

9. Definition of Terms 

di inner radius 

a 2 outer radius 

a r n , bn coefficients of the external field 

C a absorption coefficient 

C s scattering coefficient 

E magnitude of incident wave 

F frequency (Mc/s) 

ja{F, a 2 ) absorption factor in (2) 

jit* spherical functions defined by (A8) 



g relative admittance of material =-^K e /K m 

hn\ hn ] first and second spherical Hankel functions 

j n spherical Bessel function of the first kind 

k propagation constant of material 

& propagation constant of free space 

K e relative permittivity of material 

K m relative permeability of material 

W a power absorbed (watts) 

Wi power density of incident wave (watts/m 2 ) 
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W s power scattered (watts) 

z { n generalized Bessel function 

Z relative impedance of material =-y/K m /K e 

&n, In, &n dimensionless functions defined by (A4), (A5), (A6) 

e permittivity of free space =1/c 2 mo 

e* complex permittivity 

/x permeability of free space =4:ir X 10~ 7 H/m 

jjl* complex permeability 

£ s dimensionless variable for free space =k a s 

p s dimensionless variable for material =ka s 

a conductivity of material (mhos/m) 

t p unit vectors in spherical polar coordinates, p = l, 2, 3 

tan 8 e , tan 5 m electric and magnetic loss tangent, respectively. 

10. Appendix A. Derivation of the Mathematical Expressions 

Field coefficients. The following derivation is an application of the method used by Stratton 
[1941, pp. 563-569] to discuss the diffraction of a plane wave by a single sphere. The same 
solution may also be obtained by Wait's [1964] general method when specialized to this con- 
figuration. The incident plane wave in this paper has a time dependence (-\-i<at) instead "of 
Stratton 's (—iwt). The physical configuration is shown in figure 1. 

The fields in the three regions may be expressed in terms of the following spherical vector 
wave functions, where x=kr and z { n l) (x) is a given spherical Bessel function: 



m 



n, 



m _ T m M f 1 F ' n W COS * ? dP ' n Shl * 1 1 
y-*n W [± dn0 gin ^ l 2 de cog ^ 8 J, 

m __ n(n+l) (l) sin ^.lir r2 a)^i FdKsin <l>*P' n (d) cos**"] 

Zn {X)1 nW cos^ l+ xdx [XZn WJ L^cos<^ 2± ^inTsin / 3 j 



x 



In terms of these vector wave functions and the field coefficients, the fields in the three 
regions are:, 

a. Inner region (free space) 

-> OO Oft J_ 1 -> -> 

E=E e*»> 2 (_»)--^g_ [a> n m&+ib' n n?L] 

n =i n{n-f-i) 

\ Mo n=i n{n-f-l) 

where 

Z n =z Jn /t — i^O} 

b. Middle region (arbitrary material) 

E=£' e^'i:(-i)"- 7 ^^[(afma+af m ^)+{(6fn^+6fn< 3 »)] 

71 = 1 /6^-T"lJ 

ff= J^Soe'-'Z) (-9"=Sxts [(i?mffi+i. w mffi)-i(a?nffi+aM] 

\ Mo w=l ' o v '^ i~ J- y 

where 

c. Outer region (free space) 

71=1 "V& 1 i-J 
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H= 



V MO 71 = 1 



where 



n(n+l) 

^n — Jw) 6 n — ll n i* 1 — /t 0« 

The fields must satisfy the boundary conditions: 

where primes and double primes refer to t wo ad jacent regions. 

Denoting % m =k a m and p m =ka m = £ m ^K e K m , the following relations may be defined: 

»Pm ^n (Pro) Pm 

~TZ~ \_Km""n \£m)\ — Kmfon \Zm)On\<i) j On\%m) 7 (2) /> \ 7~* 
^£ra 't'ra Ism/ £ra 



(A2) 



(A3) 



(A4) 
(A5) 
(A6) 



The field coefficients, determined by substituting the fields into boundary equations, may 
be written in matrix form as: 

h«\p 2 ) 

(A7) 

Solving for a r n} (1) is obtained. The equations involving the b n 's are identical to those 
above except the relative admittance g=^K e /K m is replaced by the relative impedance Z= 

For calculational purposes the following substitutions were found to be useful: 



-^ 2) (P1) 





~o«" 







-h*>togaS>(pd 





4 2) 







Af(p 2 ) 


-W 2) (| 2 ) 


<4 3) 




— i»(&) 


&f(p 2 )S-«f(p 2 ) 


-A< 2) &)««(&)_ 


.»L 




_-j»(&)7»(&)_ 



A " ' (p) = £- / " 1) (p) ; A42) (p) = =5 j * (p) ; 



A« ) (A)=/, n (Ai)[^y.(&)-«Zi.(&)]; 
B =e «*-n>*» (pi)$ a) ( p2 ) _ e -*<P 2 -n)$(» ( ft )$w ( P2 ) ; 

r=e«'2-».»$®(p,)A»(ft)-«- i, ^'i)$»(p 1 )Af(p 2 ). 
Then the reflected field coefficient is given by 



a'„- 



h?(h) T+R 
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(A8) 



*i° (p m )=/i () Wk ffl (p m )-Zy n (U}=fn l UpJ-A l} (Pm) r^y m (U--l; (A9) 

L PmJ 



(A10) 
(All) 
(A12) 

(A13) 



For Im (p 2 — Pi)<<— 1, the second terms in (All) and (A12) are negligible compared to 
the first terms and the expression for the field coefficient becomes 

nr= j»(fe) ^n 1) (p 2 ) — Zy n (fe) rAl .N 

n Af(&)^( P2 )-Z5 w « 2 )' lAi4J 

This is referred to as the high conductivity case and for a shell thickness greater than about 1 
in. is valid for a conductivity greater than 10 4 mhos/m. 

W.= -?S?J-ib (2«H-l)(|<| 2 +|6j| 2 +Ke<+Re^) W. 
The power density of the incident wave is 

Wi=~ElJ^ W/m 2 . 

« i Mo 

Defining the scattering (or absorption) coefficient as the ratio of the power scattered (or 
absorbed) to the power density of the incident wave times the geometrical cross section, i.e., 

The dimensionless coefficients become 

O a> 

C s =f 2 S (2n+l)(K| 2 +|6j| 2 ) (A15) 

C a =-%it (27i+l)(|<| 2 +|6;| 2 +Rea;+Re 6;)- (A16) 

^2 71 = 1 

Mathematical functions. The recursion for all spherical Bessel functions, including/^, is 

, , 271—1 r f 

X 

The functions have the initial values 

/i 1) l(p)=i ^ 1} (P)=1 /i 2 l(p)=i /? > (P) = 1. 

Scattering and absorption coefficients. The scattering and absorption coefficients can be 
expressed in terms of the field coefficients by integrating the radial component of the Poynting 
vector over a sphere of radius R much larger than the outer radius a 2 . 

Separating the fields in (Al) and (A2) into incident fields, E t and H iy containing vector 
wave functions of superscript 1 and Bessel function J n , and outgoing fields, E r and H r , contain- 
ing vector wave functions of superscript 4 and Bessel function h n , the radial component of the 
Poynting vector is 
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The total power scattered, W s , and the total power absorbed, W a , are 
W a =-Re (fs* R R 2 sin 6ddd<t> 

W s =-Ke ff^ (E T eH r4 ,~E r4> H T e)R 2 sin 6dddxj>. 
Using the asymptotic forms for the spherical Bessel functions, 

j tt (x)-^sin(z-^) tf>(x)^ ^HO 
and the relationships, 

J J \sm gW sin 2 <£ sin a# cos 2 0/ 



rrc^ a co t i+^ ^ sin2 1) - ««*= 

J J \sin sin 6 sin 2 </> dd dS cos 2 </>/ 



2 7 m 2 (? ( + l) 2 , 

ior m=7i 



2^+1 
= for m^n 



one obtains 



W ' =I ir -J-S (2n+l)(|<| 2 +|6jP) W. 

#0 f Mo tt = l 
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